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Using the soft-collinear effective theory we derive the factorization theorem for the decays B — » 
M1M2 with M\ t 2 = 7r, K, p, K* , at leading order in A/Em and A/m\,. The results derived here 
apply even if a s (^inA) is not perturbative, and we prove that the physics sensitive to the EA scale 
is the same in B — > M1M2 and B — » M form factors. We argue that cc penguins could give long- 
distance effects at leading order. Decays to two transversely polarized vector mesons are discussed. 
Analyzing B — > nn we find predictions for B — > and |V u i,|/+ _ " r (0) as a function of 7. 



Decays of S mesons to two light mesons are im- 
portant for the study of CP violation in the standard 
model. In Q it was suggested that since mt,,E M 
A, nriM the amplitudes should factorize into simpler non- 
perturbative objects, and the proposed factorization the- 
orem was checked at one-loop. This approach is often 
referred to as "QCD Factorization" (QCDF). Factoriza- 
tion has also been considered in the "perturbative QCD" 
(pQCD) approach 0]. These approaches rely on a pertur- 
bative expansion in a s (EMA). The results obtained from 
factorization are quite predictive and may allow us to an- 
swer fundamental questions about the standard model. 
At the current time several important issues remain to be 
answered. These include: 1) The extent to which the re- 
sults are model independent consequences of QCD (since 
QCD is a predictive theory any model independent limit 
must give the same answer in different approaches). A 
complete proof of a factorization theorem will answer this 
question. 2) Unambiguous definitions of any nonpertur- 
bative hadronic parameters which appear are required. 
This allows the universality of parameters to be under- 
stood, as well as making clear the extent to which pre- 
dictions rely on model dependent assumptions about pa- 
rameter values. 3) Does the power expansion converge? 
If power suppressed contributions really compete with 
leading order contributions as some studies 0, Q sug- 
gest then the expansion can not be trusted. In this case 
the only hope is a systematic modification of the power 
counting to promote these effects to leading order, or an 
identification of certain observables that are free from 
this problem. 

The soft collinear effective theory (SCET) H Q pro- 
vides the necessary tools to address these issues. A first 
study of SCET factorization for B — > 7T7t has been made 
in 0. In this paper we go beyond Refs. 0,0,11] in sev- 
eral ways. We first reduce the SCET operator basis to its 
minimal form and extend it to allow for all B — > M1M2 
decays (including two vectors). Our results show that 
all of the so-called "hard spectator" contributions are al- 
ready present in the form factors, just with different hard 
Wilson coefficients. We also derive a form of the factor- 



ization theorem which does not rely on a perturbative ex- 
pansion in a s (EMA), and show that the non-perturbative 
parameters are still the same as those in the B — > M form 
factors. In our analysis long distance cc penguins lld| 
are investigated, but are left unfactorized. For natures 
values of m& and m c we give an argument why these 
contributions can be leading order. This is contrary to 
expectations that they are power suppressed 0], but in 
agreement with expectations in 0, j^, • The pres- 
ence of these contributions could introduce large LO non- 
perturbative strong phases. Even in observables that 
are free from charming penguins our results differ phe- 
nomenologically from Ref. ^ n particular while the 
power counting in Ref. requires a heirarchy in pa- 
rameters ( B7T < ( B7T , we show that SCET allows for 
other possibilities such as ~ ( Bv . We demonstrate 
that the LO SCET results are in agreement with cur- 
rent B — » 7T7T data, and find current central values favor 
Cj ^ ~ C B7r , albeit with fairly large uncertainties. 

We set M = P when discussing pseudoscalars, M = V 
for vectors, and use an M to denote either. The decays 
B -> M1M2 are mediated in full QCD by the weak AB = 
1 Hamiltonian, which for AS = reads 
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where the CKM factor is A^ /) = V pb V* f with f = d. The 
standard basis of / = d operators are (with 0\ <-> 0\ 
relative to |ll|) 



o 2 o? 



10,77,89 

i=3 



(1) 



Of 

3 ,4 
5 ,6 

7 ,8 
09,10 
077,8g 



(j>b)v-A(dp)v-A, Of = {p f -jb a )v-A(daPfj)v-A, 

{(db)v-A(qq)v-A , (dpb a )v-A(q a qp)v-A} , 

{ (db) V -A (qq ) V+A ,(dpb a ) v-A (q a 9/3 ) V+A } , 

3e — — 

-^{(db)v-A[qq)v+A , {dpb a ) V -A^ a qp)v+A}, 

3e — — 

-^(db)v-A(qq)v-A , (d/3b a ) V -A(q a qp)v-A} , 

-^Lda^{eF^gG^T a }(l+ l5 )b. (2) 



2 



b * 


c 




a 0£ s ('mvJ 


i — 





FIG. 1: Example of long distance charming penguins. The m« 
gluons are nonperturbative and LO soft gluons are exchanged 
by the b, c, c and spectator quark which is not shown. 
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FIG. 2: Example of a long distance light quark penguin which 
matches onto a power suppressed operator. The q goes in the 
n direction, the q goes in the n direction, the broken u quark 
line is soft or collinear and the u and gluon remain hard. 



Here the sum over q = u,d, s, c, b is implicit, a, (3 are color 
indices and e q are electric charges. The AS = 1 Hyy is 
obtained by replacing (/ = d) — ► (/ = s) in Eqs. 
The coefficients in Eq. are known at NLL order 
In the NDR scheme taking a s {mz) — 0.118 and m& = 
4.8 GeV gives C 7l (m b ) = -.317, C Sg (m b ) = -0.149 and 
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Ci_i (m 6 ) = {1.080 , -.177 , .011 , -.033 , .010 , -.040 , 
4.9 x 10~ 4 , 4.6 x 10" 4 , -9.8 x 10~ 3 , 1.9 x 10~ 3 } . (3) 

The relevant scales in B — ► M1M2 are m b , m c , the jet 
scale V^A and A. Varying A between 100 - 1000 MeV 
the jet scale is numerically in the range y/EA ~ 0.5 — 
1.6 GeV. Integrating out ~ mj, fluctuations, the effective 
Hamiltonian in SCETi [J^ can be written as 
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where and &^ are Wilson coefficients, the ellipses 
are higher order terms in A/Q, Q = {m;,,i?}, and Q cE 
denotes operators appearing in long distance charm ef- 
fects as in Fig. ^ Penguin contractions with light quark 
loops are included in matching onto Qff since their 
long distance contributions are power suppressed Q . The 
long-distance contributions occur when one or both of the 
quark lines in the penguin loop become soft or collinear. 
In matching onto SCET these quark lines are left uncon- 
tracted and give rise to higher dimension operators which 
are power suppressed. An example which gives rise to a 
six quark operator is given in Fig. 

In penguin contractions with charm quarks the situa- 
tion is different due to the threshold region. For the cc 
system the offshellness depends on the value of q 2 = m 2 x, 
and long distance contributions from x — > or x —* 1 are 



P 2 ~Q 2 \ 
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FIG. 3: Factorization of B -» MM' in SCET. 

suppressed U However, for q 2 ~ Am 2 the charm quarks 
are moving non-relativistically. This region corresponds 
to momentum fractions x ~ Am 2 /m 2 ~ 0.4 in the middle 
of the distribution 4>m{x). These contributions have one 
a s {2m c ), but can not be calculated perturbatively. Using 
NRQCD power counting they are "suppressed" by 0(v) 
with v ~ 0.4 — 0.5. Thus we conclude that these contri- 
butions may be leading order, and comparable in size to 
other penguin terms such as those from the small Wilson 
coefficients C3-Q. A rigorous account of these long dis- 
tance cc penguin contractions can only be obtained by 
deriving a factorization theorem for them, however we 
do not attempt to do so here, and therefore do not write 
down operators for Q cS . 

In Eq. Q the O(X ) operators are [sum over q — u,d, s] 
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with Qy^ obtained by swapping d — ► s. In Eq. (J5J the 
"quark" fields with subscripts n and n are products of 
collinear quark fields and Wilson lines with large mo- 
menta uji. For example 



n,u>2 ,^3\ 1 
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where creates a collinear quark moving along the n 
direction, or annihilates an antiquark. The b v field is the 
standard usoft HQET field with Lagrangian — b v iv- 
Db v . For a complete basis we also need operators with 
octet bilinears. We take these to be Q| 0) with T A ® T A 
color structure, for example 



? (o) = 
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^P L T A b v ] [dn, u JP L T J 
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These id and is operators do not contribute to the decays 
B — > M1M2 at leading order, but will in power correc- 
tions. Our basis of QfJ operators can be directly related 
to the one derived in except that we also included 
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Qad which makes the basis sufficient to accommodate all Here £ 



(i) 



electroweak penguin effects. 

We also need the 0(A) operators for the LO factoriza- 
tion. Defining 



they are: 
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Our basis in Eq. is simpler than the one in [£| for sev- 
eral reasons. Terms with a B^ or in the n-bilincar 
can be reduced to Eq.© by Fierz transformations. This 
shows that hard-spectator and form factor contributions 
are related. Second, fxQij — 0, so integration by parts 
allows a basis for with no n-covariant derivatives, so 
only field strengths B^ appear, plus {u n ^^PLb v ]V^[d n T^ 
PlUu) terms which give vanishing contributions. We sup- 
press Q^'s with octet bilinears that do not contribute 
at LO. The operators Qgg 1 ^ only contribute to SU(3) n 
singlet production and are not used below. 

Next we determine the most general structure of the 
p 2 ~ EA contributions in SCETi . We decouple the usoft 
modes by making the field redefinitions Q £„/ — > Y n r£ n r, 
A n ' — > Y n i A n i Y\, , with Y n i a Wilson line of n'-A us gluons 
and n' — n or n. In all Y's cancel except for 

(Y^b v ) and the operators factor into (n, v) and n 
parts, 



(10) 



In Fig. the M' meson only connects to the rest of the 
diagram at the scale p 2 ~ Q 2 , through Q™ f = qn,u 2 ^q n ,u 3 
for some flavors q, q' and Dirac structure T. The shaded 
p 2 ~ EA region is required to generate the collinear M, 
similar to the B — > M form factors At LO it is given 
by T-products of the remaining parts of the operators in 
Eq. (fTTTf) . Q^l , w hh one Lagrangian inserted on 
the spectator quark to swap it from usoft to collinear: 

Ti = fd i yd i y'T[Q\°\0),iCf] q (yUC^ n (y^ 
^C%Hyi)]+fd^T[Q?\0),*c\l: 2 q \y)], 

T 2 = fd i yT[Q ( i 1 \o),i/: ( £ q (y)]. (11) 



YigBi:W^ n + h.c. [3, and the form of 



our other C's can be found in [14j . 

Now we match SCETi onto SCET n . A complete 
treatment of Ti is an open question due to endpoint sin- 
gularities 0[ll,[ll], but (V±\Ti\B) = and the nonzero 
matrix elements can be parameterized as 



(P\T l \B)=m B C 



BP 



(V [] \T 1 \B) = m B ( 
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For T-i the most general perturbative matching at [i 2 ~ 
EA generates a set of operators with Wilson coefficients 
given by jet functions J and Jj_ whose form is constrained 
by RPI, chirality, power counting and dimensional anal- 
ysis [u>i = ZLO, UJi = (1-Z)0J, X = 1-X , \n,u = (W* 

T [(inW) Ul igB^ 4 PR,L} ia (0) [igBi-W^ n }i b (y) 

J(z,X,k + )(tf,P L ,R'y ( Lj .{Xn,xufyPL,RXn ,-xw]}, (13) 
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where {i,j} and {a, 6} are spin and color indices. At 
tree level we find that J{z,x,k + ) = J±(z,x,k + ) = 
8(x — z)Tra s (fi)CF /{N c xk + ). The remaining pieces of T 2 
are purely usoft and match directly onto soft operators 
in SCETn , giving 

-— fd A y[q s YY(y) [6^(0), (14) 
m b J 

-— fi i y[q.YY(y)[^b v ] i (Q) 1 (15) 
m b J 

where here Eq. I|14|) goes along with the J± term, and 
Eq. i|15|) goes along with the J term. 

To obtain the final result for amplitudes we combine 
Eqs. I)1()I14[I . simplify the Dirac structure between the 
soft fields, and take matrix elements. First consider 
final states containing perpendicularly polarized vector 
mesons, B — > V±V±. Kagan [l7| has argued that B — > 
Vj_Vj_ is power suppressed relative to the longitudinal po- 
larization, B -> V\\ V\\ . At LO in SCET Q { °/ 1} for i = 1-6 
have scalar bilinears and give vanishing contributions to 
B — > V±V±. The operators and generate the 
J± term in Eq. (jl3[) and could contribute. However, chi- 
rality conservation in SCETi implies that one vector is L 
and one is R polarized so the Q7/ and Qg/ contributions 
also vanish (quark masses flip chirality and in SCETi are 
suppressed by powers of m q / V AE (13). More explicitly 
the Jj_ term in Eq. I|13|l vanishes because the soft Dirac 



structure can be reduced, t^PlJ^J^ 



( a/3 



Pl, and this tensor vanishes when contracted with the 
fi-bilinear, 
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TABLE I: Combinations of Wilson coefficients appearing in the factorization formula. Note that these results do not assume 
isospin symmetry and all VV channels in this table are longitudinal. Due to our basis choice the coefficients Txj,2j(u, z) for all 
these states are identical to Tif^fC") with each c/ (tt) replaced by b- {u,z). 

Thus at LO only A cS could give transverse polarized vec- 
tor mesons so 



A(B -> V^V^ 



2Gf 
V2 



Qcc\B) . (17) 



Next consider B V\\V\\, B —* V\\P and B — > PP 
decays. Now it is the J term in Eq. I|13fl that con- 
tributes along with possible long distance charming pen- 
guins. Due to the form of our operators the J term is 
identical to the analysis of the B — ► M form factors. The 
LO factorization formula for A = (MiM2\Hw\B) which 
determines B°, B~ — ► M\M 2 with M\ >2 pseudoscalars or 
longitudinal vectors is 



A(B -> M1M2) 



\^A^ M * 



x lduT 2C (u)^ h (u) + f Ml ( 
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fB.fM 1 fM 2 



dk + J(z, x, k. 

x0 Ml (x)0 M2 (u) + Tu(u, z)^ (x)0 Ml (u) 



/m 2 C 
du T K (u) <j) Ml {u) 

) 



T2j(u,z) 
UK) 



(18) 



where A cc denote possible long distance charming pen- 
guin amplitudes which contribute in channels where c^' s ^ 
appear. For each decay mode the set of hard coefficients 
Tj£ and Tij can be obtained from Tabled 

A new result from our analysis is that the jet func- 
tion J in Eq. QlSj l is the same as that appearing in the 
factorization formula for B — ► M form factors |l9(. We 
quote here two of these formulas, one for the standard 
B — > Piv form factor f + (E), and one for the form factor 
An for B — ► Vulv decays, 



AJE): 
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msE A 2 {E) (m B + m v ) 
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where 
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2m b 



ME) , (19) 
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At LO in SCET [H E [3 E E3 
f + (E) =T(+\E)( BP (E)+N 



1 poo 

dx / dk+ 
Jo 



xC 



(+)/ 



,j ; (z, E)J(z, x, k +1 E)<f> M (x)4>+{k+) 
A ll (E) = T( A \0(E)C BV u(E) + N [l 



1 poo 

dx I dk^ 
Jo 



xCP'(z, E)J(z, x, k+,E)<t> M (x)cp%(k + ), (21) 

where N = f B fpm B /{AE 2 ),N\\ = f B fvm B /(AE 2 ), 

and the functions T(+ ,a \E),C < j~' A \z) are combina- 
tions of SCET Wilson coefficients and can be found 
in 19]. In that paper the jet functions j(- s -\z,x,k+) in 



Eq. ijT3|) are denoted by J ; 



z , x. 



and (x, 



,k + ). At the endpoint where E ~ rns/2 



the same parameters £ BM and jet function J appear in 
the form factors and in the non-leptonic decays. Since 
the analysis for J is identical to that in the form fac- 
tors several important facts can be immediately taken 
over for B — > M\M 2 decays. In particular to all orders 
in perturbation theory only the 4> B ~(k + ) wavefunction is 
obtained as proven in Ref. Also the convolution in- 
tegrals with J are finite with an identical proof to the 
one given in Ref. 0]. Finally it is clear that possible 
messenger fluctuations |2l| can not spoil factorization in 
Q,fi which have color singlet fi-bilinears, and so their 
role will be identical to that in the form factors. 

At this point we compare our result in Eq. 1|18|) with 
the result in QCDF Q. From Eq. (25) of Q the LO 
factorization theorem is 



(M 1 M 2 \O i \B) = 
{F B ^(Q)f M2 

+/mi/m 2 /b jdudxdk + 

xT} l (x, u, k + )4>Mx (x)4>m 2 {u)4>B{k+) 



(1 



(22) 
(23) 



where the 
F B ^ M (0), . 



parameters are the QCD form factors 
>Mj , and 4>b (other parameters appear when 



5 



power suppressed terms from annihilation or chirally en- 
hanced corrections are included). In the QCDF power 
counting the second term is suppressed relative to the 
first by a factor of a s . The result in Eq. ll'L'l) is quite 
similar to the SCET formula derived in Eq. (|18fl . How- 
ever, there are several important differences, which we 
comment on. The two things that are most important for 
phenomenology are that QCDF does not allow for a lead- 
ing order A™ contribution, and that the SCET analysis 
suggests that the contribution from £ and £j are com- 
parable in size, rather than (f * < ( B7r as in QCDF. As 
discussed later, current data on B — > irir seems to support 
(j* ~ ( Bv , albeit with large uncertainties. This differ- 
ence has significant phenomenological ramifications, as 
it implies that even in absence of leading order charming 
penguin effects the perturbative strong phases predicted 
in [l| would receive 0(100%) corrections. Besides these 
points there are several technical differences between the 
two formulas. Using F B ^ M (0) in Eq. rather than 
^bm d oes no t completely separate out all contributions 
from the hard scale. Also, in Eq. J23 T 1 and T n include 
perturbative contributions from both the /i 2 ~ Q 2 and 
/i 2 ~ EA scales |2(| ■ In the result in Eq. (jTB|) these scales 
are separated in T i: j and J respectively. If Q BM is inde- 
pendent of the /i 2 ~ EA scale as argued in Ref. 0] then 
the scales are also completely separated in the ( BM 
term, otherwise £ SM still encodes physics at both the jet 
scale EA and the scale A 2 . 

The jet function J depends on physics at the intermedi- 
ate scale, so its perturbative expansion in a s (V EA) is not 
as convergent as for the Tij and which are expanded 
in a s (Q). In fact perturbation theory may fail for J all 
together. This can be tested both by experiment and 
by additional perturbative calculations. Using SCET we 
can still obtain an expression for A(B — » M1M2) without 
expanding J perturbatively, 



|/ Ml jdu dz T u (u, z)C BM - (z)0 Ml («) (24) 
+fM 1 C BAh fduT K (u)^(u)\ + {l ~ 2}+\WAg M >, 



G F m 2 B 
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where power counting implies £ 



BM 



/■BM 
S,7 



(A/0) 



3/2 



Here the non-perturbative parameters ( BM , C, BM (z) 1 and 
(f> M (u) , still all occur in the B — ► M semileptonic and rare 
form factors. For a model independent analysis they need 
to be determined from data. Note that it was possible 
for us to derive Eq. I|24(l because in Eq. (|18|l we separated 
the scales Q 2 and EA into T's and J's respectively. The 
corresponding results for the form factors in Eq. I|21|) are 



/+ = T^>(E) C^{E) + N I dz C { +\z) ( BM (z, E) , 

Jo 

A,, = T^O (E) ( BV «(E) + N\\ / dzC i f ,,) (z)Cf V (z,E). 

J ° ' (25) 



The two form factors in Eq. (|21|l can be obtained from 
data on B — > (P, Vu )£v, giving important information 
on the C BM Xj M appearing in Eq. ffTHft . Note that in 
Eqs. (|T%j) and lt^4l the £'s are evaluated at E = m B /2. 
Eq. i|18|) and (|24|l are the main results of our paper. 

Using Eq. (|24|l still requires matching the full theory 
Oi's onto the Qyl' 1 ^ to determine the Wilson coefficients 



Cjf' and b\ J> . For the coefficients of Qf we find [/ = d, s] 
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V cj 3 / 27V C 

/ TO fc \ 3Cio 
V w 3 ) 2N C 

w 2 / 2N C 



(26) 



(27) 



A6<'> 



A6 



(/) 



Ab 



(/) 



-A 



(/) 
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where u>2 — rribU, 0J3 — —m^u — rrib(u — 1) and the 
Acf' and Ab\^ are perturbative corrections. The 0(a s ) 
contribution to the Ac^ (u) have been calculated in 
and later in |3| • It is possible that these results will 
need to be modified by an additional subtraction for 
the long distance charming penguin. Finally, any full 
a s {mb) analysis requires Abj (u,z) which are currently 
unknown, unless the numerical values of C, O are such 
that O - a s (m b )( so that Cf M < C BM and the Ac ( / ] 
coefficients dominate numerically. 

There are several issues in the phenomenological use 
of the factorization formula. There is a hierarchy due to 

CKM factors and the Cds which have to be accounted for 
( f ) ( f ) 

in the c\ and b\ ' . For example, C\ is about a factor of 
six larger than any of the other coefficients, making c\" , 

and large. We will refer to quantities as "con- 
taminated" if 1 /rrifc power corrections could compete with 
LO results due to the heirarchy in Wilson coefficients. 
Unless these corrections can be accounted for or proven 
to be absent, one should assign ~ 100% uncertainty to 
predictions for contaminated decays. The determination 
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of whether a quantity is contaminated depends on the rel- At tree level the b\^ 's are independent of z and this 



ative size of ( BM and (f M . If ( BM > (f M as in QCDF 
then any / — d decay in Table [i] that is independent 
of could receive large corrections, making quantities 
such as Br(B° — ► 7r°7r°) contaminated p|. Here the most 
problematic are large power corrections proportional to 
C\A/E which is ~ C 2 and 3> Ci>z. These can arise for 
example from T-products involving the operators. 

/-BM 



The situation is much better in the case ( BAI 
since any decay depending on cf"\ 6^ , or b% will not 
be contaminated and can be expected to have power cor- 
rections of normal size, ~ 20%. Our analysis of B — > hit 
below favors this situation, in which case Br(B° -» 7r°7r ) 
is not contaminated. 

At leading order in A/E there are only two sources of 
strong phases: the one- loop Aci,A&i which can become 
complex Q , and the unfactorized A c5 charming penguin. 
Additional final state phases come from power correc- 
tions ~ A/E. It is known from B° — > D°ir° decays 
that A/E corrections produce ~ 30° non-perturbative 
strong phases in agreement with dimensional analy- 
sis |22J. These large phases have nothing to do with a 
A/m c expansion so we expect strong phases of similar 
size from power corrections in B — > MiM 2 - For contami- 
nated decays, such as B — > KK, non-perturbative strong 
phases oc C\ could be order unity. 

The factorization theorems in Eqs. (|18I24|) can be 
used to make quantitative predictions for nonlcptonic 
B — » MM' decays. There are many applications; a few 
of the more important categories are: i) Decay modes 
which are independent of charming penguin contributions 
arc determined by £ and (j which can be extracted from 
semileptonic form factors, ii) SCET implies SU(3) rela- 
tions beyond those following from Hw in Eq. JTJ with full 
QCD. It also simplifies the structure of SU(3) breaking 
corrections, iii) For B — > VV SCET allows us to ana- 
lyze polarization effects, iv) Using isospin SCET makes 
predictions for matrix elements whose quantum numbers 
differ from the reduced set of A^ lM2 amplitudes. In the 
remainder of the paper we discuss examples in each of 
these categories. In particular we show that Eq. 
gives a reasonable fit to the current B — + tttt data. 

The parameters ( BM and (f M in Eq. (gljl for nonlep- 
tonic decays are common to those appearing in B — > M 
form factors Eq. (|25|) . Decays that do not depend on A c5 
include all combinations in Table [i] that are independent 
of C4 and 64, such as B~ — > 7r°7r~ and B~ —> p°p~ once 
isospin is used. For example, 

V2A(B- - n-H ) = 



(28) 



x I Jdudz (b[ d) +b { 2 d) -b { 3 d) )(u, z) (f*(z)<P*(u) 



relation gives a clean constraint on C, Bt[ and (f v = 

Flavor SU(3) symmetry is a powerful tool for study- 
ing nonleptonic B decays. In one particular application, 
Ref. proposed using flavor SU(3) symmetry to de- 
termine 7 from B + — > Ktt, tt + tt q . Corrections to this 
approach come from SU(3) breaking effects and are typ- 
ically ~ 30%. The factorization relation Eq. Ij24(l implies 
enhanced SU(3) relations beyond those in QCD. For ex- 
ample, in QCD all B — > PP decays to two pseudoscalar 
octet mesons are parameterized in the SU(3) limit by 5 
complex amplitudes. Using the SCET factorization for- 
mula Eq. Ijl 80 this number is reduced to one complex 
amplitude A c5 , one real number £ and one real function 
(j(z). In the language of Ref. [23| the operators in Eq. Q 
do not generate the E, A, and PA amplitudes, so these 
are power suppressed. 

In certain cases the SU(3) breaking can be also com- 
puted. Such an example is the determination of two 
SU(3) breaking parameters R\ t2 appearing in a SU(3) 
relation used to extract 7 |2^| 



A(B~ 



K°ir~) + V2A(B~ JWtt ) = 



(29) 



y/2^\{Ri - S EW e^R 2 )A{B- t^tt ) . 



Here Sew parameterizes the largest electroweak penguin 
effects and is calculable. The parameters Ri t2 can be 
expressed in terms of ( B7T , C, BK , (j n (z), C BK (z) and cal- 
culable Wilson coefficients and do not involve A™ or 

AKtt 
A cc ■ 

Polarization measurements in decays to two vector 
mesons have received much attention recently. These 
decays were studied in Ref. 0|j an d it was argued 
that factorization implies Rt ~ anc ^ R±/R\\ = 

1 + 0(l/mb), where Rq.t.±.\\ denote the longitudinal, 
transverse, perpendicular and parallel polarization frac- 
tions (R T = R±+ R\\, R + Rt = 1). Using SCET we 
find that Rt is power suppressed in agreement with |l7j | , 
unless the charming penguin amplitude A cS spoils this 
result. We can not resolve the validity of the R±/R\\ re- 
lation working only at LO in l/m*. Experimentally, one 
finds 0HH 



R (B+ - p+p») 
R Q (B°^p+p-) 
R (B° - cj>K*) 



= 0.975 ±0.045, 
= 0.98±g;g| zh 0.03 , 
= 0.49 ±0.06. 



(30) 



It has been argued that the large transverse polariza- 
tion observed in the (j>K* mode might provide a second 
hint at new physics in b — > sss channels beyond sin(2/3) 
from B — > 4>Ks- Unfortunately this conclusion could 
be spoiled by a contribution from A cS at leading order. 
A c5 does not contribute to B + — ► p + p°, but can affect 
B° — > 4>K* and B° — > p + p~. Until charming penguins 
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are better understood the polarization measurements do 
not provide a clean signal of physics beyond the standard 
model. 

We finally examine in some detail the predictions of 
this paper for B — > tttt decays, and show that they re- 
produce the existing data. The present world averages 
are HI] 

Sir* = -0.74 ±0.16, CW = -0.46 ±0.13, 
Br{B+ -► 7T 7r+) = (5.2 ± 0.8) x 10~ 6 , 



Br(B° -> tt+O = (4.6 ± 0.4) x 10~ 6 . 
Br(B° -> 7r°7r°) = (1.9 ± 0.5) x 10~ 6 , 



(31) 



where the branching fractions are CP averages. The am- 
plitudes are naturally divided into two pieces with differ- 
ent CKM factors, as A = \ ( u ] T + \ { c d) P, where T and 
P are usually called "tree" and "penguin" amplitudes. 
The decay amplitudes for B — > tttt can be written in a 
model-independent way as 

A(B° -» Tr+Tr-) = X^T c (l + r c e i5 "e^), 
A(B°^n\°) = AMr„(l + r „e a »^), 
\plA(B~ — > tt°tt~) = Ai d) T, (32) 

where (r c , S c ) and (r n , <5 n ) parameterize the ratio of tree 
to penguin contributions to B° — ► 7r + 7r~ and £?° — > tt q tt°, 
respectively. We have neglected small electroweak pen- 
guin contributions. Isospin gives the relations 



T = T r + T„ 



Tn^n^ 



iS n 



0. 



(33) 



leaving only 5 independent strong interaction parameters 
in Eq. 

In the first step of the analysis, we assume that (3, 7 are 
known, use this to disentangle the tree and penguin am- 
plitudes, and thus extract the five parameters in Eq. (|32ll . 
In a second step, these parameters are compared with 
the leading order predictions from SCET, and used to 
extract the nonperturbative parameters appearing in the 
factorization formula Eq. 1)24(1 . working at tree level in 
matching at the hard scale. The resulting SCET param- 
eters are then used to predict values for |V^(,|/+(0) and 
Br(B° — > 7T 7r°) as functions of 7. 

Assuming values for the CKM angles (3 and 7 we can 
use the 5 pieces of experimental data given in Eq. (|31|l to 
determine the 5 parameters in Eq. (|32|l . Using (j3, 7) = 
(23°, 64°) and the data for the CP asymmetries we find 
for the penguin parameters r c and 5 C 



= 0.75 ± 0.35 , S c = -44° ± 12° . 



(34) 



This is in good agreement with the recent determinations 
of these parameters in Refs. 27]. Using the branching 
ratio data as input, we can determine the tree parameters 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

FIG. 4: Constraints on the triangle of tree amplitudes T/T c — 
Tn/Tc = 1 from current world averaged data onB-+ rnr. The 
shaded regions show the two l-cr regions for 7 = 64° including 
the error correlation between \t\ and \t n \. The central values 
for 7 = 54° and 7 = 74° are also shown. 



as well. We find 



\T\ = AT 



(0.29 ± 0.02) ( 
1*1 = 2.07 ±0.42, \t n \ 



3.9 x 10" 



(35) 



1.15 ±0.33 (I) 
1.42 ±0.35 (II) 



where N„ = ^i r n 2 B f 7T and we defined 



T 

t = — , 



T 

t - — 



(36) 



Some of the errors in Eqs. I|34|l and l|35(l have sizeable 
correlations. The results for the tree triangle are shown 
graphically in Fig. 01 The two 7 = 64° solutions cor- 
respond to those in Eq. Q35JI and the ellipses denote la 
contours. Also shown in this figure is the isospin tree tri- 
angle, which for the reduced tree level amplitudes reads 
1 + t n = t. There are two strong phases in this triangle 
which are also shown in the figure, namely 9 between T 
and T c and 9 n between T n and T c . 

As a second step the extracted amplitudes are com- 
pared with the predictions of this paper at leading order 
in A /mi, and tree level in the SCET Wilson coefficient 
cf^ 1 and b^f\ At this order our result has four indepen- 
dent parameters. The tree amplitudes T, T c are given 
by the factorization relation Eq. I|18|l and depend on the 
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units x [ 

0.25 r 



3 .9 xW 3 ' 
\Vub\ - 




75 Y 



FIG. 5: Model independent results for ( Bw , (f n , and the 
B —* it form factor f+(q 2 — 0) as a function of 7. The shaded 
bands show the 1-a errors propagated from the B —* tttv data. 

non-perturbative parameters C B7r >Cj 7r j 



T = ^-(C 1 +C 2 )[4C^ + (4+(w- 1 )^)C^] 



(37) 



+ (d+Ci + il + iu- 1 )^ 



C-2 + C 3 



where (ir 1 )* = Jo and Cf* = /tfaCj "(4 

The penguin amplitude also gets a contribution from the 
complex A£J amplitude, so 



P 



Am 



(<*) 



c,- ( f) 



+ (C4 + (l + (^- 1 ) 7r )^)cf 7r + ^ 



(38) 



The amplitude T n is given by the isospin relation Eq. I|33|) 
&sT n = T-T c . At tree level in SCET Wilson coefficients 
the B — > 7r form factor at g 2 = is 



/+(0) = C^ + Cj 



(39) 



Neglecting the 0{a s (mb)) corrections introduces an error 
of about 10% for the T amplitudes, which is smaller than 
the expected size of the power corrections ~ 0(A/E). 

Eq. (|37|) implies that the tree amplitudes T, T c are cal- 
culable in terms of the £, £j parameters, and their rela- 
tive strong phase are small 9, 9 n ~ 0(a s (mb),A/E). On 
the other hand, the penguin amplitude P can have an 
0(1) strong phase due to the charming penguin ampli- 
tude A™. The pattern of results in Fig.QJsupports these 
predictions for the tree amplitudes T, T c for the upper 
hand solution. In particular, within the experimental 
uncertainty the phases 9 and 6 n are still consistent with 
being small and compatible with order 0(A/E) effects. 



Using the numbers in Eq. (|35|l for \T\ and \t\ and the 
SCET results in Eqs. l|3*Yjl we can extract the nonper- 
turbative parameters (,(.J- Taking LL order for the co- 
efficients (Ci = 1.107, C 2 = -0.248, C 3 = 0.011, C 4 = 
-0.025 at fi = 4.8GeV) and (u" 1 ^ = 3 0, we find 



/-Btt 



7=64° 



. (0.05 ±0.05) $*£p), (40) 
3.9 x 10~ 3 



J l 7 =64° 



(0.11 ± 0.03) ( 



where the quoted errors are propagated from the exper- 
imental errors from |r| and \t\ in En. fity. Using the 
results for r c and S c in Eq. and \V c b\ — 0.041 the 
penguin amplitude is 



P_ 



7=64° 



(0.043 ± 0.013) e 



i(136°±12°) 



(41) 



The C W7r and (T terms in Eq. JSHJl contribute 0.002 to 
P/N v , which is only a small part of the experimental re- 

f f ) 

suit. The perturbative corrections from the Aq s or 
( f) 

particularly the Ab\ s can add terms whose rough size 
is estimated to be ~ Cj^Ci ct s (mb)/n ~ 0.007. After re- 
moving these contributions, the sizeable remainder would 
be attributed to A 1 ^ . Since can have a large non- 
perturbative strong phase, the large phase in Eq. 141|) 
supports the conclusion that this term contributes a sub- 
stantial amount to P/N^. 

The extraction of the above parameters allows us to 
make two model independent predictions with only 7 and 
\V u b\ as input. First a prediction for the semileptonic 
B — > 7T form factor /+(0) is possible. Combining Eq. H40(l 
with Eq. we find 



/+(°)L=64o = (°- 17±0 - 02 ) 



3.9 x 10- 3 
I K b I 



(42) 



In Fig. we show results for ( Bn , (f* ', and / + (0) for 
other values of 7, thus generalizing the results in Eqs. i|40|) 
and Ij42(l . Note that including the correlation in the er- 
rors for £ Bv and £j T has led to a smaller uncertainty for 
/+(0). Theory uncertainty is not shown in Eq. (|42|l or 
Fig. [S] and the most important source are power correc- 
tions which we estimate to be ±0.03 on /+(0). One loop 
a s {mb) corrections are also not yet included. Varying 
fj, = 2.4-9.6 GeV in the LL coefficients changes /+(0) by 
only a small amount =f0.01. 

It is interesting to note that the central values from 
our fit to the data give £j n > £ B7r which differs from the 
hierarchy used in QCDF. Furthermore our central value 
for /+(0) is substantially smaller than the central values 
obtained from both QCD sum rules 29] (/+(0) = 0.26), 
from form factor model based fits to the semileptonic 
data |j3_f f + (0) = 0.21), or those used in the QCDF 
analysis ||j (/+(0) = 0.28 or 0.25). 
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Our analysis can also be used to make a prediction for 
Br(B° -> 7r°7r°). At tree level in SCET \t n \ = \t\ - 1 
which gives 



f (B° -> A°) 

\t\ 



(\t\-l\ 2 rj 2r c / 1 
V \t\ J |i| 2 \t\ \ \t\ 



(43) 
cos(<5 c ) cos(7) . 



Thus we predict 



Br(B° ^vrV ) 



(1.0±0.7)xl0" 6 , 
(1.3±0.6)xl0- 6 , 
(1.8±0.7)xl0~ 6 , 



7 = 54° 

7 = 64° . (44) 

7 = 74° 



These results are all in reasonable agreement with the 
current world average. The uncertainty quoted in 
Eq. (|44Jl is only from the inputs in Eq. I)43|l , and will 
be directly reduced when the first four measurements in 
Eq.(PU improve. Since the (f* term in Eq. gDJ is > 
our results for Br(B° — ► 7r°7r°) are not contaminated 
and we expect that theoretical uncertainty from power 
corrections plus a s (m&) corrections will add a ~ 20-30% 
uncertainty to the results in Eq. I|44|) . Note that one can 
turn the analysis in Eq. (|44|l around and use the data 
on B — > 7T7T in Eq. I|31(l to give a new method for deter- 
mining the value of 7, where the theoretical input from 
factorization is that the tree triangle is flat. 

Our values in Eq. I|44|l are somewhat larger than the 
central values predicted in QCDF (~ 0.3 x 10~ 6 3) or 
pQCD (~ 0.2 x 10" 6 HH). For 7 = 54° the first term 
in Eq. (|43|) dominates our result, while the r\ penguin 
term has a large cancellation with the interference term 
ex cos(7). For larger 7's this cancellation becomes less 
effective and Br(B° -> 7r°7r°) increases. In QCDF ( B7T 
dominates over a small Cj^i but has a small coefficient 
oc C2 + Ci/3, so the first term in Eq. (|43|l is small. In 
pQCD the M a ^ e terms which are multiplied by C\ are 
also small for B — ► 7r°7r°. 

In this paper we have used SCET to derive a factor- 
ization theorem for B — > M\Mi decays and explored 
the theoretical and phcnomcnological implications. Sev- 
eral issues for B — + MiM 2 still remain to be resolved. 
A factorization formula for the charming penguin con- 
tribution should be worked out, and polarization effects 
should investigated beyond leading order. It needs to 
be shown that the n—n factorization is not spoiled by 
Glauber degrees of freedom. The one loop Afro's need 
to be computed, as well as a resummation of Sudakov 
logarithms which are given by the evolution equations 
for the SCET operators. Charming penguin effects need 
to be better understood in an effective theory approach, 
and a full factorization theorem for the A cS amplitude 
should be worked out. Finally, power corrections (in- 
cluding so called chirally enhanced terms, annihilation 
contributions, and C\h./E terms) should be studied us- 
ing SCET. 



This work was supported in part by the DOE under 
DE-FG03-92ER40701, DOE-ER-40682-143, DEAC02- 
6CH03000, and the cooperative research agreement DF- 
FC02-94ER40818. IS. was also supported by a DOE OJI 
award. 



[10 

[11 

[12; 

[13; 

[14 

[IS! 
[16 
[17 

[is; 

[19 

[20 
[21 



[22 
[23 



M. Beneke, G. Buchalla, M. Neubert and C. T. Sachra- 
jda, Phys. Rev. Lett. 83, 1914 (1999); Nucl. Phys. B 591, 
313 (2000). For earlier work see Ref. 0. 

A. Szczepaniak, E. M. Henley and S. J. Brodsky, Phys. 
Lett. B 243, 287 (1990); M. J. Dugan and B. Grin- 
stein, Phys. Lett. B 255, 583 (1991); H. D. Politzer and 
M. B. Wise, Phys. Lett. B 257, 399 (1991). 

Y.Y. Keum, H.n. Li and A.I. Sanda, Phys. Lett. B 504, 
6 (2001); Phys. Rev. D 63, 054008 (2001). 
M. Beneke, G. Buchalla, M. Neubert and C. T. Sachra- 
jda, Nucl. Phys. B 606, 245 (2001). M.Beneke and 
M.Neubert, Nucl. Phys. B 675,333 (2003). 
M. Ciuchini et al., |arXiv:hep- ph/0208048 
C.W. Bauer, S. Fleming and M.E. Luke, Phys. Rev. D 
63, 014006 (2001); C.W. Bauer, S. Fleming, D. Pir- 
jol and IW. Stewart, Phys. Rev. D 63, 114020 (2001); 
C.W. Bauer and IW. Stewart, Phys. Lett. B 516, 134 
(2001); C.W. Bauer, D. Pirjol and IW. Stewart, Phys. 
Rev. D 65, 054022 (2002). 

C. W. Bauer et al, Phys. Rev. D 66, 014017 (2002). 
J.g. Chay and C. Kim. Phys. Rev. D 68. 071502 (2003); 
arXiv:hep-ph/0301262 

M. Ciuchini, E. Franco, G. Martinelli and L. Sil- 
vestrini,Nucl. Phys. B 501, 271 (1997); M. Ciuchini et 
al, Phys. Lett. B 515, 33 (2001). 

P. Colangelo, G. Nardulli, N. Paver and Riazuddin, Z. 
Phys. C 45, 575 (1990). 

G. Buchalla, A.J. Buras and M.E. Lautenbacher, Rev. 
Mod. Phys. 68, 1125 (1996). 

C.W. Bauer, D. Pirjol and I.W. Stewart, Phys. Rev. D 

67, 071502 (2003). 

M. Beneke, A. P. Chapovsky, M. Diehl and T. Feldmann, 
Nucl. Phys. B 643, 431 (2002). 

C. W. Bauer, D. Pirjol and I.W. Stewart, Phys. Rev. D 

68, 034021 (2003). 

M. Beneke and T. Feldmann , |arXiv:hep-ph/0311335| 

B. O. Lange and M. Neubert, |arXiv:hep-ph/0311345| 

A. Kagan, talks at 10 36 SLAC workshops, May & 
Oct. 2003 fc preprint UCT P-102-04; see also Y. Gross- 
man, arXiv:hep-ph/0310229 We thank A. Kagan for 
comparing results which helped us catch a sign error. 
I. Z. Rothstein, arXiv:hep-ph/0301240 A. K. Leibovich, 
Z. Ligeti and M. B. Wise, Phys. Lett. B 564, 231 (2003). 
J. W. Chen and I. W. Stewart, arXiv:hep-ph/0311285 

D. Pirjol and I.W. Stewart, Phys. Rev. D 67, 094005 
(2003); for B -> V see ibid, arXiv:hep-ph/0309053 

M. Beneke and T. Feldmann, Nucl. Phys. B 592, 3 (2001). 
T. Becher, R. Hill and M.Neubert, 
arXiv:hep-ph/0308122 T. Becher, R. J. Hill, B. O. Lange 
and M. Neubert, arXiv:hep-ph/0309227 see also 

C. W. Bauer, M. P . Dorsten and M. P. Salem, 
arXiv:hep-ph/0312302 

S. Mantry, D. Pirjol and I.W. Stewart, Phys. Rev. D 68, 
114009 (2003). 

M. Gronau et al., Phys. Rev. Lett. 73, 21 (1994); M. Neu- 
bert and J. L. Rosner, Phys. Rev. Lett. 81, 5076 (1998); 



10 



A. J. Buras and R. Fleischer, Eur. Phys. J. C 11, 93 
(1999); arXiv:hep-ph/9810260 . M. Gronau, D. Pirjol 
and T. M. Yan, Phys. Rev. D 60, 034021 (1999). 

[24] J. Zhang et al. [BELLE Collaboration], Phys. Rev. Lett. 
91, 221801 (2003); B. Aubert et al. [BABAR Collabora- 
tion], Phys. Rev. D 69, 031102 (2004). 

[25] [Belle Collaboration], Phys. Rev. Lett. 91, 201801 (2003); 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. 
Lett. 91, 171802 (2003). 

[26] H eavy Flavor Averaging Group [HFAG] : 

ht t p : / / www . slac . Stanford . edu/xorg/hfag/ 
[27] A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, 



Phys. Rev. Lett. 92, 101804 (2004); C. W. Chi- 
ang, M. Gronau, J. L. Rosner and D. A. Suprun, 
|arXiv:hep-ph/0404073| A. Ali E. L unghi and 
A. Y. Parkhomenko, arXiv:hep-ph/0403275 

[28] A. P. Bakulev, S. V Mikhailov and JM. G. Stefanis, 
|arXiv:he p-ph/0303039 

[29] P. Ball and R. Zwicky, JHEP 0110. 019 (2001). 

[30] Z. Luo and J. L. Rosner, Phys. Rev. D 68, 074010 (2003). 

[31] C. D. Lu, K. Ukai and M. Z. Yang, Phys. Rev. D 63, 
074009 (2001). 



